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1 The Colored Glass Condensate 

At very small values of the Bjorken x-variable, one expects QCD to enter a new regime which 
is caracterized by parton saturation and very high values of the QCD field strength i* 1 " ~ 1 jg. 
Saturation, which is a limitation on the maximum phase-space parton density that can be 
reached in the hadron wavefunction, may have been already observed at HERA0, and should 
play an important role in the very early stages of relativistic heavy ion collisions at RHIC and 
LHCu (and Refs. therein). In the saturation regime, the individual parton-parton interactions 
may be weak |] (which we shall assume in what follows; i.e., we assume that g <C 1), but the 
parton densities are so large that the system becomes strongly non-perturbative. Thus, at 
a theoretical level, understanding saturation is a challenging and fascinating problem where 
one has to deal with fully non-linear QCD. This is reminiscent of a similar problem in high 
temperature QCD where perturbation theory breaksJown at the soft scale g 2 T because of the 
large thermal occupation numbers of the soft gluons Q (and Refs. therein). 

The efforts toward understanding the region of high gluon density have uncovered a new 
form of matter which is formed from these gluons u. This matter is universal in that it should 
describe the high gluon density part of any hadron and nuclear wavefunction at small x. The 
combination of high density and the fact that gluons are massless bosons leads naturally to the 
expectation that this matter is a Bose condensate. Since the gluons carry color and local color is 
a gauge dependent quantity, any gauge invariant formulation will neccessarily involve an average 
over all colors to restore the invariance. This averaging procedure bears a formal ressemblance 
to the averaging over background fields done for spin glasses 0. The matter is therefore called 
the Colored Glass Condensate (CGC). 

It should be emphasized that this picture holds, strictly speaking, only in the infinite- 
momentum frame, where the hadron propagates almost at the speed of light, and thus appears 
as an infinitesimally thin two-dimensional sheet (by Lorentz contraction). In this frame, the 
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Figure 1: A parton cascade leading to a small- a; virtual gluon. 

parton interpretation makes sense and deep inelastic scattering (DIS) proceeds via the instan- 
taneous absorbtion of the external probe (e.g., a virtual photon 7* with 4-momentum q^) by 
some parton in the hadron. The Bjorken xb parameter is defined as xb = Q 2 /2P ■ q, where 
Q 2 = —q^q^i and P^ = 5 fJ/+ P + , with large P + , is the hadron 4-momentum fl. By kinematics, xb 
coincides with the longitudinal momentum fraction x = p + /P + of the struck parton: xb = x. 

At x <C 1, the gluon density increases faster, and is the driving force toward saturation 
|]. The dynamics that leads to this increase is the quantum evolution toward small- x : in a 
parton cascade initiated by some fast parton (i.e., a hadron constituent with a relatively large 
longitudinal momentum k + ~ P + , e.g., a valence quark), the number of radiated gluons increases 
exponentially with the rapidity gap At = \n(k + /p + ) ~ ln(l/x) between the original parton and 
the soft (i.e., small p + : p + = xP + <C P + ) final gluon which is struck by the external current (see 
Fig. 1). This (BFKL) picture — which assumes the radiated gluons to behave as free partons — 
ceases to be valid at very small x, where the gluon density is so large that the radiated gluons 
overlap each other in the transverse plane. This is the onset of saturation. 

This is also the regime where the description in terms of a colored class condensate becomes 
appropriate: Corresponding to the strong ordering in longitudinal momenta along the cascade, 

k + = k£ > kf » k^ » • • • > k^ = p + , (1) 

there is a similar hierarchy among the lifetimes of the radiated gluons (since the latter are 
proportional to the respective k + momenta): Softer a gluon is, shorter is its lifetime. 

This hierarchy in time stays at the basis of a quenched approximation Iffl! for the "fast" 
degrees of freedom: when "seen" by the soft modes with very short lifetimes, the modes with 
larger p + appear as frozen (no dynamics) and can be replaced by a static (i.e., independent of 
x + ) and random color charge configuration, with density p a (x~,x±). (This is random since the 
soft gluons can belong to different cascades.) The spatial correlations of the effective charge 
Pa(x~,x±) reflect the quantum dynamics at large longitudinal momenta, and are encoded in a 
statistical weight function W[p]. Because of the Lorentz contraction, we can write p a (x _ ,xj_) 
5(x~)p a (x±), and W is a functional of the superficial charge density p a (x±) alone. 

Thus, in order to compute soft correlations in this approximation, one has to first study 
the (quantum) dynamics of the soft gluons in the presence of a given color charge p a , and then 
perform a (classical) average over p a , with weight function W[p]. Clearly, the latter will depend 
upon what we call "soft" and "fast", i.e., upon the separation scale A between fast (p + > A) 
and soft (p + < A) degrees of freedom. We thus write W[p] = Wa[p\- 

c We use light-cone vector notations: for some arbitrary vector p M , we write p M = (p + ,p~ ,p±), with p + = 
(l/V2)(p°+p 3 ), p~ = (l/V2)(p°-p 3 ), and p ± = (p\p 2 ). 

d To directly measure the gluon density, it is convenient to consider a Gedanken experiment where the DIS is 
initiated by the "current" j = —\F^"F^ U which couples directly to gluons. 



For instance, the 2-point correlation function is obtained as (in the light-cone gauge A + = 0) 

with the functional integral running only over soft gluon fields A^(p) withp + < A, and the action 
S[A, p] describing the dynamics of these fields in the presence of the classical source p a , in the 
eikonal appoximation 0: S = Sym + Sw, where Sym = J d 4 x(— F^ v /A) is the usual Yang-Mills 
action, and (with x = (x~,x±), and T denoting time-ordering of color matrices): 

Sw = jfjr J d 3 xTr |/9(x)Texp ig J dx + A~(x + ,x) j . (3) 

The double averaging in eq. (||) is very similar to the one performed for spin systems in a random 
external magnetic field (which plays the role of p in the above equation), and supports the 
physical picture of the saturation regime as a coloured glass condensate. 

Of course, the final results for soft correlators must be independent of the arbitrary sepa- 
ration scale A. The A-dependence of the weight function Wa[/°] must cancel against the cutoff 
dependence of the quantum theory for the soft modes. This constraint can be formulated as a 
renormalization group equation for W\ [/)]□, to be presented in Sect. 3 below. 

2 Saturation in the classical approximation 

Consider first the simple approximation where the quantum path integral in eq. (|2|) is evaluated 
in the saddle-point (or classical) approximation SS/dA^ = 0, and the weight function is taken 
simply as a Gaussian (this is the McLerran-Venugopalan model Q, originally formulated for a 
large nucleus for which the Gaussian approximation is expected to work better) : 



w A [ P ] - ex P j-^r / dV^(z±)j. 



(4) 



Here, p\ is to the total color charge squared (per unit area) of the partons with p + > A. 
The classical approximation of eq. @/ reads (with x = {x~,x±)) : 



(A^x + ,x)A»(x + ,y)) d = J V P W A [p]A»(x)A»(y), (5) 

where A^(x) is the solution to the classical Yang-Mills equations with source p a (%), 

[D V ,F^] = gS^S(x-) Pa (x ± ), (6) 

and is time- independent, like p a itself. The LC gauge solution can be written as0: 

A + = A- = 0, A\x) = e{x~)A\{ Xl _) + 0(-x-)AL(x±), (7) 

with A l + (x±) and A l _(x±) related to p(x±) via a non-linear equation. Note that the vector 
potential A l (x~ ,x±) is discontinuous at x~ = 0, so the associated electric field F +% = d + A l is 
localized at the light-cone (i.e., within the support of the source): 

T +i {x) = S(x-)(A^(x ± ) - AUx±)). (8) 

By using this approximation, let us compute the gluon distribution function, that is, the 
number of gluons per unit of x in the hadron wavefunction having transverse momentum less 
than Q. This is defined as (with = (A; + ,kj_) and k + = xP + ) : 

xG(x,Q 2 ) = 1 1 |^e(Q 2 -& ± )(ltf(z+,k)lF(x+,-k)). (9) 



In the classical approximation F l+ (x + ,k) ~ J- l+ (k±) is independent of both x + and k + , 
and (with A A 1 = A\ — AL, and the classical average defined as in eq. @) : 

xG d (x,Q 2 ) = ^jcfx ± e- ik ^(AA i a (0 1 _)AA i a (x ± )) d . (10) 

Here, A is the hadron transverse aria, and we have assumed transverse homogeneity. Thus, 
xG c i(x, Q 2 ) is independent of x, which reflects the absence of quantum evolution in the present, 
classical approximation. With the Gaussian weight function (|4|), and the non-linear classical 
solution (0), the classical gluon distribution (|To| ) can be computed exactly One obtains: 

, > AT 2 - 1 1 - e - x ± ln (^Q C D)Q 2 s/^ 

(A4i(0J AA^)) d = ? , (11) 

where Q s oc a s fi\ is the saturation momentum, and is a priori a function of A. Remarkably, this 
equation displays saturation via non-linear effects in the classical solution. This interpretation 
can be made sharper by going to momentum space. If N(k±) is the Fourier transform of (pi]), 

N(k±) oc a s (Q 2 s /k 2 ± ) for kj_ » Q 2 , (12) 

which is the normal perturbative behaviour, but 

N{k ± ) oc — In % for k\ « Q 2 S , (13) 
a s Qi 

which shows a much slower increase, i.e., saturation, at low momenta (with kj_ ^> Aqcd though). 
According to eq. (|Tl]), saturation is also a statement about the maximum field intensity that 
can be reached in the system: the classical field never becomes larger than A 1 ~ 1/g. This is the 
maximal occupation number for a classical field, since larger occupation numbers are blocked 
by repulsive interactions of the gluon field. 



3 The non-linear evolution equation 

In writing down the effective theory for soft gluons in eq. @ , we have assumed that the influence 
of the fast gluons can be reproduced by a classical color source p with a peculiar structure (time- 
independent, and localized at x~ = 0), and some (still unspecified) weight function Wa[p]. In 
this section, we show how to construct this effective theory, step by step, by integrating quantum 
fluctuations in successive layers of p + . 

To this aim, it is convenient to consider a sequence of two effective theories ("Theory I" 
and "Theory II") defined as in eq. (0), but with different separation scales: A in the case of 
Theory I, and bA for Theory II, with b-C 1. That is, Theory II differs from Theory I in that the 
"semi-fast" fields with longitudinal momenta between bA and A have been integrated out, and 
the associated correlations have been incorporated at tree-level, within the new weight function 
WbA- The difference AW = Wb\ — Wa can be obtained 111 by matching calculations of soft 
{k + <C bA) gluon correlations in both theories. The result can be expressed as an evolution 
equation for WV[/°] (with r = ln(P + /A)) with respect to variations in r. 

The quantum corrections due to the semi-fast fields have to be computed to leading logarith- 
mic accuracy (LLA), that is, to leading order in a s ln(l/6) — indeed, it is only to this accuracy 
that the hierarchy of scales in eq. ([l]) is satisfied, and the matching is possible — , but to all 
orders in the classical fields and sources (since, in the saturation regime of interest, the non- 
linear effects are so strong that cannot be expanded in perturbation theory). This specifies the 
accuracy of the evolution equations to be obtained. 



To this accuracy, there are only two contributions to AW : the one- and two-point correlators 
of the fluctuating colour charge 5p a (x) of the semi-fast gluons. Specifically, one obtains illl 

(5p a (x)) p a s \og{l/b)8{x~)G a (xA_), 
(5p a (x)5p b (y)) p « a s log(l/6) 5(x~)xab(x±,y±) 5(y~), (14) 

while all the n-point correlators with n > 3 are of higher order in a s . In these equations, (•••)/> 
denotes a quantum expectation value over the semi-fast fields at fixed p. Furthermore, a a (x±) 
and Xab( x ±yU±) are generally non-linear functionals of p(x±) given by one-loop diagrams within 
Theory I (with loop momenta restricted to the strip: 6 A < \p + \ < A). 
In terms of these functions, the evolution equation for W T [p] reads "fli 

^4^*'' -^4 (i5) 

This functional equation is equivalent to an infinite hierarchy of ordinary equations for the 
correlators of the charge density. For instance, by multiplying eq. ([0|) with p x Py and functionally 
integrating over p, one obtains an evolution equation for the two-point function: 

"J- (PxPy)r = a s (Xxy + Px^y + V X py)r , (16) 

ur 

which in general, however, involves also the higher n-point functions, via a and x- But in the 
weak source approximation, i.e., with a and x computed to lowest order in p, this becomes 
a closed equation for (pp) which has been shown Q to be equivalent to the BFKL equation, as 
necessary on physical grounds. This is a highly non-trivial check of the effective theory in eq. (|2|). 

In order to study saturation, however, one needs eq. (|lj|) in the regime of strong background 
fields and sources (A 1 ~ 1/g and p ~ 1/g 2 ; cf. eqs. ( pT| ) and @), which requires for an exact 
calculation of the coefficients a and x- This has been done recently a, via a lenthy calculation 
which had to cope with difficulties related to gauge-fixing, the axial poles in the gluon propagator, 
and the proper definition of the singular limit p(x) —* 5(x~)p(x±). This makes possible to look 
for solutions to eq. ([l5|). Note that, formally, this is like a functional Schrodinger equation in 
imaginary "time" r. An interesting possibility is that the "Hamiltonian" in its r.h.s. has an 
eigenfunction W[p] of maximum eigenvalue A. This would lead to an universal behavior of W r [p] 
in the small-x, or high-energy, limit: W r [p] — > T -*oo e TX W[p], with r-independent W[p]. 
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